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ABSTRACT 

We describe a search for Mg II AA2796, 2803 absorption lines in Sloan Digital Sky Survey (SDSS) spectra of 
QSOs whose lines of sight pass within impact parameters p ~ 200 kpc of galaxies with photometric redshifts 
of z = 0.46-0.6 and errors Az ~ 0.05. The galaxies selected have the same colors and luminosities as the 
Luminous Red Galaxy (LRG) population previously selected from the SDSS. A search for Mg II lines within a 
redshift interval of ±0.1 of a galaxy's photometric redshift shows that absorption by these galaxies is rare: the 
covering fraction is f(p) ~ 10-15% between p = 20- 100 kpc, for Mg II lines with rest equivalent widths of 
W r > 0.6 A, falling to zero at larger p. There is no evidence that W r correlates with impact parameter or galaxy 
luminosity. Our results are consistent with existing scenarios in which cool Mg Il-absorbing clouds may be 
absent near LRGs because of the environment of the galaxies: if LRGs reside in high-mass groups and clusters, 
either their halos are too hot to retain or accrete cool gas, or the galaxies themselves — which have passively 
evolving old stellar populations — do not produce the rates of star formation and outflows of gas necessary to 
fill their halos with Mg II absorbing clouds. In the rarer cases where Mg II is detected, however, the origin of 
the absorption is less clear. Absorption may arise from the little cool gas able to reach into cluster halos from 
the intergalactic medium, or from the few star-forming and/or AGN-like LRGs that are known to exist. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: halos — galaxies: groups: general — 
galaxies: clusters: general — quasars absorption lines 



1. INTRODUCTION 

Red galaxies contain th e bulk of all the stel lar mass that ex- 
ists in the local universe ( Hogg et al. 2002 ), and thus mark 
the conclusion of both the dominant star formation processes, 
and of the assembly of galaxies. The Luminous Red Galax- 
ies (LRGs) are a population selected by their optical colors 
to contain old, passively evolving stellar populations at inter- 
mediate redshift. The first two principal samples of LRGs 
were cons tructed from the Sloan Digital Sky Surv ey (SPSS; 
York|2000| ) for galaxies at redshifts z ~ 0.2-0.4 (Eisenstein 
et aLpOOlj, and from the 2d F-SDSS LRG And QSO (2SLAQ) 
survey ( [Cannon et al.|20 06 ) for galaxies atz^0.4-0.7. Sub- 
sequently, the properties of the LRGs have been investigated 
by many authors: LRGs hav e luminosities L > L * and large 



stellar masses M> 10 n M o (Ban erji et al.|2010| ); the major- 
ity of LRGs show spectral energy distributi ons (SEDs) with 



little evidence for on-going star formation (Roseb oom et al. 
2006); they resi de in dark matter (DM) h alos with masses 
M ~ 10 13 - 14 M« JMandelbaum et al.||2006|), and in Halo Oc- 



cupation Distribution (HOD) models, most LRGs are ce ntral 
galaxies dominating the DM halos (e.g. Zheng et al. 2009] and 
refs therein); they appear to be most prevalent in clusters (Ho 
|eraLl[2009] ). At z ~ 0.2-0.8, the evolution of the LRG lu- 
minosity function (LF) is found to be consistent with a mode 
in which stars form at high-z, and e volve passively thereafter 
with little further star formation (e.g. | Wake et a l. 2006; Brown] 
|et al.|2007| |Cool et al |20Q8l |Banerji et al |2010j ). It appears 
that these massive galaxies have already been assembled by 
z ~ 1 ( [Conselice et al.|2007|[Ferreras et al. 2009; Banerj i et al.| 

Given the colloquial view of LRGs as "red and dead", they 
would seem an unlikely population of galaxies to harbor sig- 
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nificant amounts of cool gas. Nevertheless, there has recently 
been considerable interest in using the clustering of LRGs to 

ption lines arise 
Gaut hier et al.| 



infer the masses of ha l os in which Mg II absor pi 
dBouche et al.||2006t [Tinker & Chenj[2008l [( 



2009; Lundgr enet al.||2009| — see jp~2| for a detailed dis 



cussion of these results). The Mg II AA2796,2803 doublet is 
the first strong UV transition to be redshifted into the opti- 
cal wavelength regime, at z > 0.2 — a redshift low enough 
for galaxies responsible for the absorption to be identified. 
Historically, initial searches for galaxies at the same redshifts 
as previously i dentified Mg II systems were q uick to confirm 



as previously identified Mg 11 systems were quick to confirm 
an association (|Bergeron|[l9 86; Yan ny et al.| 19871 [Cri stiani 
T987I |Bergeron & Boisse|199lV|Le Brun et al.J_1993 ). More 



extensive surveys suggested that all galaxies with luminosities 
~ L* have absorbing cross-sectio ns of radii R s ~ 60 kpojwith 
unity cove ring fraction f c (e.g. |Steidel|[T993| |1995[ |Steidel| 



et al.|1994| ). The subsequent development m numerical hydro- 
dynamical simulations, and the refinement of semi- analytical 
modelling, has provided a natural explanation for multi-phase 
gas in galactic halos, through the inflow of baryons into re- 
gions of high density. Particularly appealing for the origin 
of Mg II systems are those models which describe cold flows 
along DM filaments, where unshocked gas at temperatures of 
< 10 5 K are channeled into galaxies (e.g. |Keres et al.| 2009 ; 
|Brooks et al.||2009] |Kaufmann et al.||2009| and refs therein) 
at the same time that shocked gas is accreted by the halo. As 
Mg II absorption is thought to arise in photoionized gas with 
temperatures ~ 10 4 K, the existence of halos with Mg II ab- 
sorbing clouds is consistent with these simulations. However, 
one of the more important parameters governing the existence 
of cool gas in DM halo models is the mass of the halo. The 
simulations predict that in halos with masses > 1O 13_14 M , 

3 Throughout this paper, we adopt Hq = 70 km s _1 Mpc -1 , Q m = 0.3 and 
= 0.7. Distances cited from published works that are originally given in 
units of h~ x kpc are converted to this cosmology. 
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where the LRGs reside, there should be no cool gas, whereas 
in lower mass halos, where the L* galaxies mentioned above 
are found, there could still be enough cool gas to cause Mg I I 
absorption along sightlines that intercept the halos (see §4.2j . 
In these scenarios, QSO sightlines passing close to LRGs 
should show no Mg II absorption. 

On the other hand, recent studies have focused less on the 
relationship of Mg II absorption and halo mass, and more 
on the idea that Mg II systems originate in some combina- 
tion of galactic disks and starburst driven outflows. Galac- 
tic outflows offer a possible (although observationally un- 
substantiated) explanation for finding cool gas clouds out 
to many tens of kpc from a galaxy. The fact that large- 
scale outflows are driven by significant star formation activ- 
ity has led to more detailed investigations int o the link be- 
tween Mg II systems and late-type galaxies. |Bouche et al.| 
(2007), for example, showed that two-thirds of 21 sampled 
Mg II systems with rest-frame equivalent widths (REWs) of 
W r > 2 A have Ha emission from galaxies with star forma- 



tion rates of 1 -20M© yr _i . Zibetti et al. (2007) found ex 



cess (broad-band) emission around stacked SDSS QSO im- 
ages with known Mg II systems, compared to stacked images 
of QSOs without absorption. The optical colors of the light 
suggested that the SEDs of strong Mg II absorbers are bluer 
than for weak ones, indicative perhaps of the domination by 
blue star- forming galaxies to the production of strong lines. 
Moreover, the large sky coverage of the SDSS, combined with 
the remarkable number of QSO spectra obtained, has made 
it possible to examine the relationship between star forming 
galaxies and Mg II systems, without regard to the identifica- 
tion of any galaxy close to the line of sight. Emission lines 
have been detected in QSO spectra at the same redshifts as 



absorption line systems, ei ther in stacks of spectra (Wild et al. 
[2007; Menard et al. 20 10]) or in rarer cases, individua l QSO 
spectra ( |Noterdaeme et al.|2010t jBorthakur et al |2010 



Whether these results indicate that absorption arises in 
large-scale gaseous outflows, rather than simply from galac- 
tic disks that host star forming regions, is less clear. The idea 
that galactic dis ks are responsible for Mg II lines has long 
been appeali ng (Wagon er|1967HBowen et al.| 1995 [Cha rlton 



& Churchill||1996HSteidei et ai.||2002| ), and modeling Mg II 

"isk 



absorbing gas as a combination of disk and outflow features, 
naturally explains the range of galaxy luminosities an d star 
formation rates observed (Chelouche & Bowen 2010). The 
detection of low ionization optical absorption lines in out- 
flows from nearby starburst g alaxies, using the light from the 

rriiliiYA/ itc<=»1f to niv^hp fh<=» craQ flVTiirtin 1 QQQ ■ ^r'Viw'irtv Rr 



galaxy itself to probe the gas ( Martin 1999; Schwartz & Mar 
tin|2004t|Martin|2005j[Rupke & Veilleux|2005| ) adds support 



for the model, and indeed, Mg II absorptio n is seen from out- 
flows of high-z galaxies (Tremo nti et al.||2007| |Weiner et al.| 
|2009| ). Detailed investigations of starburst galaxies certainly 
sugge st that these types of galax ies produce multi-phase out- 
flows ( [Strickland et al.|2004bja| ) that would show absorption 
lines from many different ionization states if penetrated by a 
QSO sightline. It remains to be seen though how ubiquitous 
these outflows are at higher redshift around the majority of 
late-type galaxies, whether QSOs lines of sight, probing in 
a direction transverse to galaxy sightlines, would see similar 
absorption signatures as those recorded toward the centers of 
galaxies, or, indeed, whether outflowing gas capable of pro- 
ducing the types of absorption seen can escape out to many 
tens of kpc from galaxies. 
The redshifts at which Mg II becomes detectable in SDSS 



spectra [and where the SDSS sensitivity allows for detection 
of weaker lines in higher signal-to-noise (S/N) data] are in the 
range z ~ 0.36-0.43 (3800-4000 A). An interest in the as- 
sociation between Mg II systems and LRGs arises from the 
fact that at these redshifts and beyond, the LRGs are, prag- 
matically, the only galaxies detectable in SDSS images with 
reliable photometric redshifts (hereafter "photo-zs"). The 
availability of good quality photo-zs can be used to cull a 
large number of faint, intermediate-z galaxies without hav- 
ing to first measure their redshifts through expensive, time- 
consuming spectroscopic surveys. If the galaxies are observed 
as part of a survey — such as SDSS — which also records the 
spectra of a similarly large number of QSOs, then the opportu- 
nity to study the absorbing properties of the galaxy ensemble 
is obvious. The fact that, in this particular case, such a set of 
galaxies happens to be LRGs, means we can explore Mg II 
absorbing gas in the halos of a restricted, specific, and rare 
type of galaxy and their environment. The theoretical con- 
siderations discussed above suggest that we should find no 
absorption, but this paper seeks to test that prediction. 

Most surveys designed to find galaxies responsible for ab- 
sorption systems have been based on targeting Mg II systems 
known ab initio. While this approach yields information on 
the origin of the absorption, it does not necessarily define 
the absorbing properties of galaxies. Prior to the launch of 
the Hubble Space Telescope (HST), efforts to study the gas 
in disks and halos of galaxies focused on pairing very low- 
redshift galaxies with QSO sightlines and searching for ab- 
sorption lines (Ca II AA3933,3968, Na I AA5 889,5895, etc.) 
that could be observed from the ground (e.g. Boksenberg & 



Sarg ent] [1978, 
Bergeron et al. 



B oksenberg et al.||19801 |Blades et al.| 



1981 
991). 



1987[|Wombleet al. 1990; Bow enetaLf 
With HST, it became possible to search for the more sensitive 
and abundant UV lines from nearby galaxies instead. That 
work showed that at z ~ 0, there is little evidence for Mg II 
absorbing clouds around galaxies beyond ~ 50 kpc ( |Bowen| 
|et al.||1995| ), a nd when detected, the absorption originates in 
galactic disks (Bah call et al |1992| [Bowen et al.|1995[ [T996] ) 
or in tidal debris b etween the Ide ntified galaxy and an often 
fainter companion ([Bowen et al.||1994] |1995|) — even when 
the targeted galaxy is a known starburst i |Norman et al.|1996] ). 

Unfortunately, studies of galaxies in the local universe are 
hampered by the necessity of having to use QSOs with the 
brightest UV continua, in order to obtain sufficient S/N with 
the available HST spectrographs. To avoid this limitation, and 
to make use of the huge number of QSO sightlines observed 
spectroscopically with SDSS, more aggressive attempts have 



been made recently [following a much earlier study by Bech 
|told & Elling son] ( |1992| )] to identify galaxies at higher red 



shifts and to again search for Mg II lines shifted into the op- 
tical, without any prior information on the absorption lines in 
the background QSO spectra. The absorbing cross-section of 
a galaxy, R g , and the covering fraction of the gas f c , are ob- 
viously inter-related; if f c declines with radius (i.e. if there is 
no sharp boundary to R g ), then R g can be set arbitrarily for a 
given f c , and vice versa. Similarly, both R g and f c can, in prin- 
ciple, depend on the equivalent widths of the lines. In a study 
using early results from SDSS, Tripp & Bowen (2005 ) found 
galaxies at z ~ 0.3-0.6 with f c only ~ 50% for R g ~ 50 kpc, 
for sensitivities of W r > 0.1 A. Similarly, for W r > 0.3 A 
ton &Cooke| ([2009]) found f c 
I 



at z 



1TI5 
0.1. 



0.25- 

for smaller radii of R„ 



Bar- 
kpc. 



0.4 for R g ~ 110 
■ g ^ 50 kpc, for gal axies 
Conversely, however, Chen et al. ( 2010| ) [see 
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also Chen & Tinker ( 2008| )] measured much higher cover- 
ing fractions and larger cross-sections for galaxies at z — 0.2: 
R g ~ 110 kpc and f c -[0.7,0.8] for W r > [0.3,0.1] A. These 
radii are obviously considerably larger than the cross-sections 
measured for galaxies at z ~ 0. 

One caveat to this work is that at higher redshifts, identify- 
ing the processes responsible for absorption systems is much 
harder. At these redshifts, galaxy features are hard to discern 
even in high resolution HST images, low surface brightness 
features are difficult to detect, multi- wavelength observations 
(e.g. 21 cm interferometry observations, which are important 
for mapping high H I column densities around galaxies) are 
often not possible, and the light of a background QSO masks 
any faint emission directly along a sightline. (Indeed, these 
difficulties also ham per the study described in this paper, as 
we discuss in §4.1| ) Hence, even with more comprehensive 
surveys of absorbing galaxies at z > 0.1, it remains difficult 
to know whether Mg II absorption marks gas inflowing from 
filaments, outflows from star formation activity, the extension 
of galactic disks, the tidal debris between multiple galaxies, 
or some combination of all these mechanisms. Hence our 
current ambitions are to pre- select galaxies with defined char- 
acteristics, whose properties have been determined from un- 
derstood processes, and use such calibrations to interpret the 
origin of the observed absorption systems. 

In this paper, we study the absorbing properties of LRGs 
— again, selected without any prior bias toward knowing 
whether or not Mg II exists at the redshifts of the galaxies. 
Our initial approach was to simply define a sample of SDSS 
galaxies close to QSO sightlines that had photometric red- 
shifts with very small errors, and to search for Mg II absorp- 
tion in the SDSS spectra of the background QSOs. As de- 
scribed below, we began such a program without regard to the 
type of galaxies selected. We found, however, that the selec- 
tion of LRGs in the SDSS arises naturally from demanding 
that the photo-z errors of galaxies be small. We note here 



that, concurrent with our investigations, Gauth ier et al.| |2010 ) 
have also examined the absorbing properties of a set of galax- 
ies (with photo-zs) selected a priori to be LRGs. Their results, 
obtained independently of ours, obviously offers an important 
comparison with our study, and we discuss their results in $3] 
and fj] 

The outline of this paper is as follows: in j p.l| we define 
how we selected the galaxies to study, and de scrib e th e sea rch 
for Mg II at their photometric redshifts in |2.2| In |2.3| we 
describe the properties of the galaxies, and demonstrate why 
the majority of them are LRGs. We discuss the results of the 
study in ^3] and their implications in ^4] 

2. SAMPLE SELECTION 

2.1. Selection of galaxies with photometric redshifts 

We began an investigation of galaxies close to SDSS QSOs 
with no regard to the properties of the galaxies themselves. 
That is, there was no initial aim of selecting only LRGs in 
our study. We aimed to utilize both the photometric redshift 
information that was available for SDSS galaxies, and exist- 
ing high quality spectra of SDSS QSOs. Only after forming 
a suitable sample of QSO-galaxy pairs did it become clear 
that the galaxies were LRGs. We disc uss t he properties of the 
galaxies in our sample more fully in |2.3| In this section, we 
describe the criteria that were actually used to select galaxies 
close to QSOs. 

We first compiled a list from DR6 of QSOs observed spec- 



troscopically by SDSS that had galaxies within 30 arcsec of 
the QSO sightlin^] and that also had cataloged photometric 
redshifts. We imposed a basic cutoff for the photo-z of a 
galaxy to be Zcc2> 0.4 (see below for a definition of this quan- 
tity), Zcc2<£qso> and a galaxy morphology ("type" = 3). The 
minimum photo-z of Zcc2= 0.4 was imposed since Mg II was 
not expected to be readily detectable below ~ 3900 A in SDSS 
spectra (see below). This basic selection yielded a sample of 
-80,000 pairs. 

The most important criterion in selecting a suitable subset 
of galaxies for which we would search for Mg II absorption, 
was that the photo-z should be reliable — that is, the photo-z 
error should be small. The SDSS database contains t wo mea- 
sures of a ga laxy's photo-z estimated from the work of Oyaizu 
eTaT] ( [2008] ), tagged in the SDSS database as "Dl" and "CC2 
(which we refer to as zdi and Zcc2), both neural network based 
estimators, zdi uses galaxy magnitudes in the photo-z fit, 
while zcc2 uses only galaxy colors. Each of these photo-z 
estimates have errors associated with them (based on an em- 
pirical "Nearest Neighbor Error" method), although there are 
two measurements of the error in each case, because the dis- 
tribution of errors Az = |z P hot _ £spec | — the difference between 
the photo-z and the real redshift of a galaxy in a training set 
of galaxies — are non-Gaussian. There are two estimates in 
the error of the photo-zs: one error is given by a^, which 
contains 68% of the values of Az in the training sample; the 
other, which we call cr rms (simply "<r" in Oyaizu et al.), is the 
rms scatter (again, see Oyaizu et al. for more details). 

The accuracies of CC2 and Dl decline quickly for r-band 
galaxy magnitudes fainter than r = 20.0. The majority of 
galaxies found close to QSOs in our initial search for QSO- 
galaxy pairs had r > 20.0, and excluding these would result in 
a sample too small to be of interest. Based purely on galaxy 
magnitude, the errors in zdi and zcc2 rise sharply from 0.05 
to 0.09 between 20.0 < r < 21.0 (see Fig. 5 of Oyaizu et al.). 
Fortunately, for the galaxy redshifts we are interested in — 
at z > 0.4 — photo-zs are defined to a high accuracy as the 
4000 A break in the galaxy spectra passes through the r-band. 
Figure 6 of Oyaizu et al. shows that between z = 0.4-0.6, a^s 
for zcc2 ranges from ~ 0.04-0.09, and is smaller than a rms 
by Az « 0.02. In this paper, we adopt as the error in a 
galaxy's photo-z, and denote this simply as a(z). 

We interpret Oyaizu et al.'s results to mean that, for 68% of 
selected QSO-galaxy pairs, the true redshift of the galaxy lies 
between Zphot-^fe) and z p hot + crfe)- From our original colla- 
tion of galaxy pairs, we therefore selected only galaxies with 
cr(z) < 0.05. We took the maximum photo-z to be 0.6, while 
the minimum photo-z was not determined by the photomet- 
ric properties of the galaxies themselves, but by the minimum 
wavelength so that Mg II could be detected in SDSS spec- 
tra. We took this minimum to be 3800 A, plus a redshift path 
which allowed a search for Mg II within 2<j(z) (see below). 
The resulting minimum redshift was 0.46. Given this redshift 
range, and the 30 arcsec limit for finding galaxies, the maxi- 
mum impact parameters studied are between 175 and 201 kpc. 

To further exclude galaxies with spurious photo-zs, we also 
demanded that |zcc2~^di | < o-(z), i.e. that the photo-zs esti- 
mated from the two different methods agreed to within la of 
each other. Additional constraints included: Zqso - £cc2> 0.3, 
to avoid Mg II associated with the QSO itself; and zqso< 2.1, 

4 That is, objects which are identified as "neighbors" to QSOs in the 
specObj view of the SDSS database 
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Figure 1. Representative normalized spectra of QSOs in our sample with 
sensitivities W\{ m < 0.4 A. Spectra are selected for QSOs that have the 14 
closest galaxies, with their impact parameters p printed at the top right of 
each panel. The x-axis shows the wavelength reduced to a redshift range for 
Mg II A2796, and normalized to zcci- The QSO designation is given at the 
top left of each panel. Mg II, when detected, is marked with two vertical 
bars. The dotted lines at Az = ±0.1 represent the range over which Mg II 
was searched for. 

to exclude high-z Lya-forest lines in the SDSS spectra. In a 
few cases, more than one galaxy satisfying all these criteria 
was found in a single QSO field; these galaxies were flagged 
and are discussed later. 

Further culling of the original sample was made based on 
the spectroscopic properties of the background QSOs. These 
are discussed in the next section. 

2.2. The search for Mg II absorption from foreground 
galaxies 

For the QSO-galaxy sample described above, formed from 
criteria applied to galaxies, we next considered a set of cri- 



teria based on the suitability of the QSO spectra. We first 
rejected all QSOs with S/Ns in the continuum of less than 6.0 
per pixel at zcc2, based on the error arrays generated from the 
pipeline extraction of the one-dimensional SDSS spectra. At 
z ^ 0.5, this S/N corresponds to a 3a REW limit of W r « 0.7 A 
(when calculated using the method described below). We re- 
jected spectra with extraction defects, Broad Absorption Line 
complexes, complicated continua, or those contaminated by 
complex and numerous higher-z systems. In some cases, we 
rejected spectra where Mg II was expected to fall close to a 
particularly sharp QSO emission line, since defining the con- 
tinuum at the peak of these emission lines can be difficult (par- 
ticularly in cases where associated absorption lines are seen 
superimposed on the emission line). 

For the accepted spectra, we normalized the QSO continua 
by defining a lea st- squares fit of Legendre polynomials to the 
continuum flux ( Sembach & Savage |1992| ). For each spec- 
trum, we then calculated an array of observed 2a (Wi) EWs, 
where, at each pixel i, a 2 (Wf) = YTi ( a i ^) 2 > an d is the 
wavelength dispersion in A pixel -1 . Error arrays from the 
SDSS pipeline again provided values of a t , and we took n 
to be a value twice the resolution expected for SDSS spectra^] 
or five pixels. 

In each normalized QSO spectrum, we searched for Mg II 
absorption within a window of ±2a(z) of zcc2 , i.e. Az = 
±0.1. Using a value of ±2a(z) corresponds conveniently to 
assuming that w 95 % of the LRG photo-zs will be correct; 
if we had chosen to use a^ instead of cr 68 , we would have 
to have increased the path length by Sz ~ ±0.04 to include 
±2a rms . Reviewing the SDSS QSO spectra used in our sur- 
vey indicates that one more Mg II system would have been 
detected if we had adopted this estimate of the photometric 
redshift errors. As we shall see below, this would have made 
little difference to the final results. 

Detection of Mg II required that both members of the 
doublet be found with EWs > 2a(Wj), although only the 
Mg II A2796 line had to lie in the ±2cr(z) window. When 
detected, we fitted Gaussian profiles to the lines and calcu- 
lated a final redshift for the absorption system by averaging 
the redshifts measured for each line of the doublet. Line EWs 
were measured in the standard way, W = YTi (1 wnere 
Ft is the normalized flux at the i'th pixel, and the sum is cal- 
culated over n pixels. As noted above, for unresolved lines, 
n = 5. However, many strong Mg II lines found in SDSS spec- 
tra have widths that are larger than this; in order not to under- 
estimate W therefore, for lines detected with EWs > 6a(Wi) 
(i.e. lines defined with high S/N), we increased n to cover a 
width of ±3cr g auss, where a g2LUSS was the width of the Gaussian 
profile fitted to a line. 

Were the redshift of a galaxy to be measured spectroscop- 
ically, the 2a EW limit for Mg II if a doublet was not de- 
tected would be given by the value of 2a (Wi) at the appropri- 
ate wavelength. In order to define an EW limit over a wave- 
length range corresponding to ±2cr(z), however, we need to 
be sensitive to lines with EWs greater than some EW limit 
over the entire range of ±2<r(z). We therefore define the 2a 
REW limit, Wu m , to be the maximum value of 2a (Wi) in the 
EW array between ±2a(z), after each point in the EW array 
is corrected to be in the rest frame of Zi = Aj/Ao - 1 (where Ao 
is the wavelength of the Mg II A2796 line). For SDSS spec- 

5 At the wavelengths where Mg II is expected the resolution of the SDSS 
spectra is ~ 2.2 A at 4500 A, FWHM, or ~ 150 km s _1 . 
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tra, a(Wi) is relatively flat as a function of wavelength over 
the redshift range considered here, so is close to the value that 
would be calculated at exactly zcc2- 

Examples of the spectra used in our survey are shown in 
Figured For this figure, the QSOs chosen are sorted by QSO- 
galaxy impact parameter, and we have only selected spectra 
with sensitivities of Wu m < 0.4 A, simply to highlight better 
quality data. Figure[2]shows color representations of the fields 
around these QSOs. Each image is a cutout extracted from the 
SDSS database, and so is a g-, r-, and /-ban d composite of 
the original imaging data (Lupton et al. 2004 ). Full details on 
all the pairs selected are given in Table [lj the table includes 
details of the QSOs and galaxies in the sample, as well as the 
results from the search for Mg II absorption: values of REWs 
when Mg II was detected, or limits when the lines were not 
detected, are listed therein. 

2.3. Galaxy Properties 

In Figure [3] we plot the g- r and r—i colors of the galaxies 
in our sample. Even though the colors of galaxies were not 
used in their selection, their distribution in the color diagram 
is confined to a narrow region with 1.5 < g—r < 2. and 0.7 < 
r—i < 1.2. T hese colors are simi l ar to t hose used by Eisenstein 
|eTaLl ( |200T1 ) and |Cannon et al.| ( [2006] ) to define LRGs at the 
same redshifts as the galaxies in our sample. The reason that 
these authors applied such color cuts initially, and the reason 
for expecting their surveys to focus on z ~ 0.2-0.5 early type 
galaxies, was because the SEDs of such galaxies are predicted 
to fall above specific color cuts. So, for example, Figure [3] 
shows the cuts, designated d±_ and cm, used by Cannon et aL 
to define spectroscopic targets in the 2SLAQ sample. They 
defined d± to select early type galaxies at increasingly higher 
redshifts, and cy to eliminate late-type galaxies altogether. In 
the figure, we have drawn a line for d± =0.55, which defined 
"sample 3" of the galaxies studied by Cannon et al. and which 
contained the bulk of their LRG sample. Most of the galaxies 
we selected lie in the same color region as their LRGs. 

The way in which the SEDs of early-type galaxies produce 
different broad-band colors as a function of redshift is well 
understood. A change in the g — r color to a constant value 
at g — r ~ 1.7 measures the transition of the 4000 A break in 
a galaxy's SED as it moves from the g-band to the r-band 
at z ~ 0.4; g-r begins to increase again at z = 0.7 as the 
4000 A break moves into the /-band. To demonstrate this 
effect, we hav e used the LRG SED provided by Blan ton &| 
Roweis ( 2007 ) to show how the colors change with z when 
convolved with SDSS filters. These are shown as red circles 
in Figure [3] where a circle is plotted every & = 0.02 up to 
z = 0.6. We have also plotted the expected c olors of a late- 
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type Sbc and Scd galaxy using the SEDs of |Coleman et al. 



( 1980), as green and orange circles, respectively. 

Figure [3] shows that the colors of the galaxies in our sam- 
ple are largely consistent with early-type LRG colors, and not 
with late-type galaxies. More sophisticated modeling shows 
that theoretical SEDs can also reproduce these colors, assum- 
ing stars formed in a single burst 1 1 Gyr (z ~ 2.5) ago and 
with a star formation timescale of 1 - 3 Gyr (Banerji et al. 
2010, and refs. therein). 



Hence, the majority of the galaxies in our sample are early- 
type LRGs. The reason for principally selecting LRGs is al- 
most certainly because these types of galaxies, with a well 
defined 4000 A break, provide the most accurate photo-z es- 
timates. Since we selected galaxies based on small values 
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Figure 2. Representations of some of the SDSS QSO-galaxy pairs used 
in this survey. The fields shown include those with QSO spectra plotted in 
Figure [T] they thus represent the pairs with the smallest impact parameters, 
and with Wu m < 0.4 A for the QSO spectra. Each image is 60 arcsec on a 
side, with NE to the top-left. Although not labeled, the impact parameters 
between QSO sightlines and galaxies are shown in Figure 1 and are between 
20.0 (top left corner image) and 40.0 kpc (bottom right corner image). Full 
details can be found in Table [T] The QSO designation is printed at the top 
of each image, and its position is indicated by cross-hairs. Cross-hairs in red 
represent QSOs that show Mg II near the photo-z of the galaxy, while white 
cross-hairs indicate that no absorption was detected. Galaxies are identified 
by arrows, and always lie at the center of the field. 



of a(z), we therefore predominantly picked galaxies with the 
strongest B aimer breaks. 

We can compare the luminosities of the galaxies in the 
sample with those of the LRG population. By adopting the 
same LRG SEDs discussed above, we can derive suitable K- 
corrections to apply to our galaxies and thereby estimate accu- 
rate absolute magnitudes M. These corrections are necessary, 
since the values needed to correct, for example, M r derived 
from r-band observations, and the true rest- frame M r , are 
K r ~ 0.9-1.4 mags for z = 0.46-0.6. The luminosity func- 
tion ( LF) for LRGs a t intermediate redshift has been derived 
by | Wake et al.| ( |2006] ) and is reproduced in Figure]?] The blue 
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Figure 3. Plot of broad-band colors for galaxies in our selected sample of 
QSO-galaxy pairs. Black circles represent galaxies for which Mg II was de- 
tected within 5z±0. 1 of the photo-z of the galaxy, zcci\ g ra Y circles represent 
non-detections of Mg II. To avoid confusion, errors in the colors are shown 
only for Mg II detections, and not for Mg II non-detections. Also shown are 
theoretical colors for a Luminous Red Galaxy (LRG, in red), an Sbc (green) 
galaxy, and an Scd (orange) galaxy, calculated over redshift intervals of 0.02. 
The last points plotted are for z = 0.6, and for emphasis, the redshifts corre- 
sponding to the criterion used to select galaxies, zcc2= 0.46-0.6, also have 
black outlines around the circles. The color cuts used by Cannon et al. ( 2006 1 
are shown as dotted lines, and are cm = 0.7(g-r) + 1.2(r-i"-0.18) > 1.6, and 
dj_=(r-i)-(g-r)/8.0 > 0.55 (see text). 

points represent the volume density of LRGs at0.5<z<0.6 
taken from the 2SLAQ, corrected for passive evolution as well 
as including ^-corrections. In this case, the corrections to the 
absolute magnitudes Mo.2r are made to measure the luminos- 
ity at z = 0.2. Making the same correction, but not including 
an evolutionary correction (which is much smaller than K r ), 
we measure the distribution of Mo.2r for the galaxies in our 
sample and plot them as a histogram in Figure [4] The actual 
normalization of the histogram is not relevant here: the figure 
simply demonstrates that the galaxies have largely the same 
distribution of luminosities as the LRG population. 

For the final sample of galaxies, we searched for existing 
(DR7) SDSS spectra of the galaxies, in order to confirm both 
their photometric redshifts and the shape of their SEDs. Only 
two spectra were available, for two galaxies at distances of 
more than 144 kpc from their background QSO sightlines, 
neither of which were associated with any Mg II absorption 
lines. The S/N of both spectra were very poor, and insufficient 
to provide accurate redshifts. Given the limit of the SDSS 
MAIN galaxy sample, r = 17.8, and the extended LRG sam- 



ple, r= 19.2 (Eisenstein et al. 2001), we would expect very 



few of our galaxies to have been observed spectroscopically: 
our sample consists of only [0,0.3]% of galaxies brighter than 
r =[17.8,19.2]. 

3. RESULTS 

Having applied all the selection criteria discussed in §2.1| 
the final sample consisted of 593 QSO-galaxy pairs. There 
are 26 galaxies that have another galaxy with a photo-z lying 
closer to the QSO sightline, or, 23 galaxies that have another 
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Figure 4. Plot o f the r-band l umino sity function of LRGs normalized to 
z = 0.2 taken from Wake et al. ( 2006) (blue line). In order to compare the 
absolute magnitudes of the galaxies in our QSO-galaxy pairs sample, we 
overplot a histogram (light gray) line for all the galaxies selected. The nor- 
malization of the histogram is arbitrary, but shows that the galaxies in our 
sample have similar luminosities to the LRG samples. The distribution of 
absolute magnitudes of galaxies that can be associated with Mg II absorption 
are plotted in dark gra y. The top x-axis s hows the ratio of L/L* assuming 
M* +51og/z = -20.4 (Blanton et al. 2003) and no evolution in the LRG LF 
between z = 0.1 -0.2. 

galaxy with a photo-z lying closer to the QSO sightline and 
for which the background QSO spectrum has W\ im < 0.6 A. 
For the remainder of this section, we exclude galaxies fur- 
ther away than the identified galaxy (with a pho to-z) , but 
we reconsider if their exclusion is warranted in |4.1| Of 
the 567 remaining galaxies, [175,42] have impact parameters 
of p < [100,50] kpc. At these separations, we found Mg II 
absorption within Az = ±0.1 of the galaxies' zcc2 in [41,9] 
cases. Representative spectra have already been shown in Fig- 
ure [T] Of the 89 detections of Mg II absorption (at all radii), 
11 spectra showed 2 or 3 Mg II systems within the search 
window. This number is too small to infer anything about the 
clustering of Mg II systems in this sample, and more compre- 
hensive discussions on the cross-correlation between Mg II 
syste ms and LRGs have b een given by |Bouche et al.| ([2006 ) 
and Tinker & C hen ( 2008). In this paper, we select the system 
with the redshift closest to the photo-z of the LRG as the one 
"associated" with the galaxy or its environment. 

The REWs, W r , for detected Mg II A2796 lines, or 2a 
REW limits Wn m when no lines were detected, are plotted 
against galaxy impact parameter p in Figure [5] For each p, 
there is an uncertainty a(p) caused by not exactly knowing 
the redshift of the galaxy. This error is small though, since 
\a(p)\ = 2pa(z)/(l+z) 2 , which is only 7% of p for a(z) = 0.05 
at z = 0.5. These errors are omitted in Figure [5] The fig- 
ure shows that the EW limit of the survey is < 0.7 A, as ex- 
pected from selecting QSO spectra with a minimum S/N ra- 
tio. That is, we were able to always detect any Mg II lines 
with W r > 0.7 A. In many cases of course, the spectra were 
sensitive enough to detect lines much weaker than this, down 
to W r > 0.2 A. Also as expected, galaxy impact parameters 
extended out to « 200 kpc, the limit imposed by selecting 
galaxies within 30 arcsec of a QSO sightline. Mg II lines are 
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Figure 5. Mg II A2796 eest equivalent width W r , or 2cr limits, plotted against 
QSO-galaxy impact parameter p. 
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Figure 6. Plot of the covering fraction f(p) for galaxies in the sample, plot- 
ted against QSO-galaxy impact parameter p. The gray squares represent f(p) 
without any corrections for the chance detection of an Mg II system along 
a sightline; the fraction of detections to total number of pairs in each bin 
[which defines f(p)] is written next to each point. The black squares repre- 
sent f(p) after subtracting the predicted number of chance detections, assum- 
ing dN/dz = 0.3atz~0.5. Horizontal error bars show the adopted bin sizes. 
If some of the absorption detected does not actually arise from the probed 
LRGs, but, e.g., from companion galaxies, then the points draw n here (and 
in Fig -IT} represent only upper limits to f(p) for LRGs (see g4.1) . Inset: The 
ratio oftne observed number of Mg II systems to the expected number, again 
assuming dN/dz = 0.3. 

detected up to REWs of W r ~ 2 A, but there is no correlation 
of W r with p. In particular, there are many sightlines that show 
no Mg II to REW limits of 0.2-0.7 A. Figure [^demonstrates 
that Mg II detections are not a function of galaxy color, since 
detections are found for all values of g - r and r - i colors in 
the sample. 

The covering fraction of galaxies, f(p), is defined simply 
as the number of galaxies with detected Mg II at some im- 
pact parameter, compared to the total number of galaxies. To 
calculate f(p), it is necessary to use an EW limited sample 
of sightlines. To that end, we selected only sightlines with 
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Figure 7. Plot of the covering fraction f(p) for different values of Wn m , 
after correcting for the predicted number of Mg II absorbers along the QSO 
sightlines. For the Wn m = 0.3 A sample alone (shown with black points), the 
number of detections to the total number of sightlines selected is indicated 
as a fraction next to each point. The legend in the top right corner matches 
the color of the points to the W\[ m adopted in calculating f(p), and shows 
the values of dN/dz used to correct the observed fraction of detections to 
non-detections. 

S/Ns high enough to reach REW limits < 0.6 A, irrespective 
of whether Mg II was detected or not. This meant that the 
selected spectra could detect lines as strong as, or stronger 
than, 0.6 A. Since the highest S/N spectra in such a sample 
could also detect lines weaker than 0.6 A, we only counted 
a detection if it too had a REW > 0.6 A. Hence the cover- 
ing fractions include sightlines where we could have detected 
lines with REW > 0.6 A, and where W r > 0.6 A when a 
line was detected. The results are shown in Figure [6j where 
values for f(p) are plotted as gray squares. Error bars are 
taken to be simply Poissonian errors on the number of de- 
tections. The covering fraction is consistently low, starting 
at 0.2 between p = 24-50 kpc, and declining to 0.1 between 
p= 150-200 kpc. 

Although we detect Mg II within ±2a(z) of zcc2 towards 
some galaxies, in order to properly measure f(p), we must 
also consider how many of the Mg II detections arise simply 
from chance. The numb er of Mg II absorbe rs per unit red- 
shift, dN/dz, is given by|Nestor et al.|(|2005]), and at z = 0.5 



is ~ 0.3 for W r > 0.6 A. From their work, we estimate an er- 
ror in dN/dz ~ 0.05, and include this in computing the error 
bars in Figure [6] Over a given redshift interval, we can there- 
fore calculate tne number of Mg II absorbers expected to arise 
by chance and subtract these from the number of detections. 
The resulting values are shown as black squares in Figure [6] 
The ratio of the observed number of systems to the expected 
number is also shown inset to Figure [6] The figure shows that 
f(p) is small but non-zero at p < 100 kpc, but that beyond 
this, f(p) is consistent with zero. This implies that any de- 
tection of Mg II beyond p > 100 kpc is probably just from 
chance, and that the galaxies selected in this sample have no 
associated Mg II absorption beyond 100 kpc. 

We initially selected an REW limit of 0.6 A for two rea- 
sons: first, because the limit lies close to the median of the 
REW distribution in our sample, it would allow us to include 
a significant number of spectra; and second, because a value 
of dN/dz for that limit was available in the literature. In 
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their sear ch for Mg II from a sa mple of galaxies defined to 
be LRGs, Gauthi eFet al.| p010) calculated a covering frac- 
tion for lines with a REW of 0.5 A, a value that is similar 
to ours. They found f(p) ~ 20-40% for p < 100 kpc, and 
a largely constant (with some deviations) f(p) of ~ 20% be- 
tween p ~ 100-400 kpc. This result is clearly different from 
our measurement of a sharp decline of f(p) to approximately 

z ero beyond 100 kp c. 

|Gauthier et al.| ( |2010| ) also calculated f(p) for much 
stronger lines, those with REW limits > 1.0 A, and found a 
factor of w 2 smaller f(p) than for their analysis using REWs 
limits of 0.5 A. We can see whether any differences in f(p) 
exist for different REW limits in our sample. Since we are pri- 
marily interested in f(p) after correcting for Mg II lines that 
might arise by chance in our Az window, we need to know 
dN/dz at these REW limits. For W\j m greater than 0.6 A , we 
can use values of dN/d z computed by Nestor et al. (20051 and 
|LundgrenetaL] ( [2009| ) of dN/dz = 0.2 and 0.18 lor W^T= 0.8 
and 1.0 A, respectively. Again, f(p) depends weakly, but non- 
negligibly, on the error in dN/dz, so we have estimated an 
uncertainty of « 0.05 in dN/dz. The results are shown in Fig- 
ure [7j For Wii m = 0.6-1.0 A, f(p) is largely the same, reach- 
ing between 10 - 15% for p < 100 kpc, then falling sharply at 
larger p. 

We can also estimate f(p) for Wn m < 0.6 A. Nestor et al. 
measure dN/dz — 0.7 for Wn m = 0.3 A, although the error on 
this value appears to be larger than for the other values of 
dN/dz — we estimate an error of « 0.1. Figure [7] shows that 
at this REW limit, f(p) may be perhaps as much as a fac- 
tor of two larger at p < 100 kpc, although it too falls to zero 
quickly at p > 100 kpc. Unfortunately, we have far fewer 
spectra in our sample with a S/N that satisfies the limit of 
W\[ m = 0.3 A. The errors are correspondingly large, and we 
cannot be sure that the larger f(p) values arise merely as a 
consequence of poor statistics. To emphasize this point, we 
label in Figure [7 the actual fraction of Mg II detections to the 
total number sightlines in each bin (as also shown in Fig. [6]) 
for this particular Wu m . So for < p < 50 kpc, there are only 
13 sightlines available, and we detect Mg II in five of them. 
In the interval 50 < p < 100 kpc, we find Mg II in 11 of 33 
sightlines. In both cases, the errors are large and compatible 
with the same values of f(p) found for higher Wu m . 

Searching for Mg II lines near a REW limit of 0.3 A intro- 
duces potential biases. At high REW limits, the Mg II dou- 
blet ratio, DR= (W r X2196)/(W r A2803), of detected systems 
is expected to remain low, close to unity, since both lines are 
likely to be (close to being) saturated. This means that if a 
A2796 line is detected at some significance (in our case, at a 
2a level), then the A2803 will often be detected at a similar 
significance. At a low W\{ m of 0.3 A, however, the Mg II DR 
is likely to be higher (anywhere between 1.0 and 2.0), and a 
A2796 line with a strength close to Wn m may be detected, but 
the A2803 line may not In our selection process, we re- 
quired that both members of the doublet be detected at a 2a 
significance level, in order to reject lines from other species 
at different redshifts. This selection process could, in princi- 
ple, lead us to under-estimate the number of weak systems at 

6 Mg II DRs are actually difficult to accurately measure in these SDSS 
data, since the fractional errors on both of the W r lines are relatively high, 
leading to large errors in DR. Although there is a clear trend for systems with 
W r > 0.6 A to have DRs closer to 1 .0, the number of detections for W r > 0.3 A 
is too low to provide a good comparison of DRs between the two samples. 
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Figure 8. Plot of Mg II A2796 against the absolute magnitude M r of the 
galaxies in our sample. Galaxies at impact parameters of p < 100 kpc are 
plotted in red, while the remainder are shown in black. AT-corrections were 
calculated for each galaxy using the LRG SED template discussed in %23\ 
(assuming no evolution) and zcc2- 

low REW limits with W r ~ Wn m . This would then increase the 
difference between f(p) for Wn m = 0.3 A and for higher Wn m 
searches in Figure [7] 

A comparison or f(p) for different REW limits with those 
found by Gau thier et al.|([2010| ) is complicated by the fact that 
the REW limits may be defined slightly differently for the dif- 
ferent analyses; for example, the number of pixels over which 
a line is measured may be different for each analysis. Al- 
though there is a suggestion in our data that f(p) changes by 
a factor of 2 as Wn m changes from 0.3 to 0.6 A, there is no 
change in f(p) for any increase in W\[ m when Wn m > 0.6 A. 
Gauthier et al also see an increase in f(p) by a factor of 2, but 
for different REW limits, 0.5 and 1.0 A. It is possible that we 
are seeing the same change, and that the values of Wu m are dif- 
ferent simply because of the differences in the way REWs are 
measured. We find it hard to understand, however, why the 
differences in Wu m should be so large, particular since both 
analyses use SDSS error arrays and 2a limits. We also con- 
tinue to see the sharp decline in f(p) beyond 200 kpc, which 
Gauthier et al do not see for both Wn m = 0.5 and 1.0 A. 

Figure [8] plots the REW of detected Mg II absorption lines, 
or limits, against the absolute magnitude of the galaxies. In 
this case, the absolute magnitude is ^-corrected to a redshift 
of ze ro (as opposed to a rest frame of z = 0.2 discussed in 
§2.3| ). The figure shows no clear relationship between the lu- 
minosity of a galaxy and the REW of Mg II when detected. 
The figure also shows the Mg II detections segregated by p, 
with galaxies associated with Mg II detections at p < 100 kpc 
shown as red circles, and galaxies with associated absorption 
at p > 100 kpc shown as black circles. Again, there is no 
obvious difference between the two groups. 

The distribution of the difference between the photometric 
redshift of the galaxy, zcc2, and the absorption redshift z a bs of 
the Mg II system, Av = cAz/(l+Zcc2), is shown in Figure [9] 
We do not expect this distribution to have a high enough ve- 
locity resolution to distinguish between absorption from indi- 
vidual LRG halos, or absorption from lower mass halos that 
cluster around the larger mass LRG halos (for example). In- 
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Figure 9. Histogram of the difference between the photometric redshift of 
the galaxy, zcc2> and the absorption redshift z a bs of the Mg II system, when 
detected. The gray histogram shows the distribution for all sightlines where 
Mg II is detected, whereas the orange histogram shows the distribution for 
those pairs with p < 100 kpc. 

deed, although there is some suggestion that the values of Av 
tend to cluster around zero km s -1 , any peak to the distribution 
is poorly defined, as we would expect given the photo-z error 
associated with zcc2- A plot of Av against p (not shown here) 
also shows no convincing evidence that Av increases with p. 

4. DISCUSSION 

In this paper, we have constructed a sample of QSO- 
galaxy pairs from SDSS DR6 based on galaxies having pho- 
tometric redshifts of between z = 0.46-0.6, and with small 
photometric-redshift errors of ~ 0.05. These photo-zs are 
taken directly from the SDSS datab ase, and are based on work 
described by Qyai zu et al.| ( |2008| ). Due largely to selecting 
galaxies with small photo-z errors, as well as needing to use 
galaxies with magnitudes near the detection limit of the SDSS 
imaging data, we have selected a sample made up primar- 
ily of LRGs, as defined by their colors and luminosities. We 
have searched for Mg II absorption in the SDSS spectra of the 
background QSOs over a window Az = ±0. 1 of the photo-z of 
the galaxies. Our most significant result is that the detection 
of a Mg II system by an LRG is rare, even at small impact 
parameters of p < 100 kpc. A simple calculation of the cov- 
ering fraction as a function of p suggest that the probability of 
detecting strong (W r > 0.6 A) lines within « 20-100 kpcQof 
an LRG is only 10-15%. 

4. 1 . The Role of Unseen Interlopers 

Interpreting the covering fraction f(p) requires some care. 
The most obvious question is whether or not the covering frac- 
tions shown in Figures 6 and 7 represent only lower limits to 
the true values of f(p). For a galaxy at any particular p, if 
Mg II is not detected, then the result is clear: the absence of 
absorption is completely valid for the LRG, even if there ex- 
ists another galaxy at a similar redshift which lies closer to 
the QSO sightline. Under the assumption that the photo-z of 
the galaxy is correct, then the absence of absorption from that 
galaxy is unambiguous. 

7 This lower limit of 20 kpc is set by the SDSS's ability to separate QSOs 
and galaxies in the imaging data; 20 kpc at z = 0.5 correspond to an angular 
separation of 3.3 arcsec. We implicitly assume that LRGs have not been 
probed at radii less than 20 kpc in this work. 



However, when Mg II is detected, our selection process 
does not rule out the possibility that another galaxy lies closer 
to the sightline and is actually responsible for the absorption. 
Conventionally, surveys to find absorbing galaxies are con- 
ducted with the aim of identifying all galaxies in a field, to 
some stated magnitude limit and impact parameter. In our 
survey, selecting galaxies based solely on having small errors 
in their photo-zs does not take into account the possibility that 
a galaxy exists closer to the sightline, and may be causing the 
absorption. Late-type sub-L* galaxies are not easily detected 
in SDSS imaging data at the redshifts of the LRGs: assum- 
ing M r * = -21.2 ( [Blanton et al^|2003] ) a late-type L* with a 
flat continuum below a rest-frame wavelength of 4000 A (for 
which a negligible ^-correction would be required) would 
have an observed magnitude ofr~21atz = 0.5, about 1 
mag brighter than the detection limit of the SDSS. A visual 
inspection of the fields around the QSO-LRG pairs in the 
SDSS data shows that the LRGs are often accompanied by 
very faint galaxies with intermediate colors, some of which 
may even lie closer to the QSO than the LRG. Unfortunately, 
these objects have either no photo-zs, or photo-zs with very 
large errors, making it difficult to generalize about their asso- 
ciation with the identified LRGs. In some fields, the faintest 
objects are not identified in the SDSS database, and even 
brighter galaxies are not detected if they overlap the QSO it- 
self. An example of this can be seen in Figure |2} although the 
LRG Jl 12013.16+671 116.2 lies only 38 kpclrom the QSO 
Jl 12012. 11+671 115.9, a galaxy can be seen due south of the 
sightline which is overlapping the QSO sightline. The galaxy 
is not cataloged by SDSS, and no photometric information is 
available. Mg II is detected at the redshift of the LRG, but 
whether this, or the interloping galaxy is responsible for the 
absorption, is not known. The galaxy may be unrelated to 
the LRG, or could be interacting with it. Without more pre- 
cise redshift information for all the galaxies in the field, and 
deeper, higher resolution images, an understanding of this ab- 
sorption system is far from complete. Finally, we note that 
the presence of other faint galaxies close to a given QSO- 
LRG pair is not simply confined to fields where Mg II is 
detected. Fields where Mg II is not detected are sometimes 
equally complicated. 

If other, gas-rich galaxies closer to the QSO sightlines than 
the LRGs are indeed responsible for detected Mg II absorp- 
tion, then our association of absorption systems with the LRG 
would, obviously, be erroneous and would provide incorrect 
information on the absorbing characteristics of LRGs. In Fig- 
ures]^ and [7] therefore, we would need to remove such sight- 
lines from the appropriate p bin. As this would lead to a lower 
fraction of detected absorption to total sightlines, we would 
infer that the observed f(p) were only upper limits to the true 
value of f(p) for LRGs. 

Our inability to identify potential interlopers is the neces- 
sary trade-off needed to study large numbers of galaxies us- 
ing their photometric redshifts instead of completing expen- 
sive redshift surveys of many individual fields. Instead of 
considering which galaxies are nearest to a sightline, an al- 
ternative approach to analyzing our data might be to simply 
estimate the probability of finding Mg II absorption around 
any LRG, regardless of its environment. The plots of f(p) in 
Figures!^ and [7] effectively measure this probability, if we add 
back allthe LRGs that were excluded because another LRG 
had been identified closer to the sightline (^3]). In fact, leav- 
ing in these LRGs makes practically no difference to f(p); of 
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the 23 LRGs that have an LRG closer to a line of sight (and 
for which Wu m < 0.6 A) only five show associated Mg II, and 
their impact parameters are between 90-183 kpc, a range of 
p that is dominated by non-detections of Mg II. This small 
number of detections adds a negligible amount to f(p) over 
this range. 

4.2. Possible Reasons for the Lack ofMg II Absorption 

As discussed in §1, LRGs are believed to inhabit some of 
the highest over-dense regions in the universe. In Halo Oc- 
cupation Distribution (HOD) models, they reside in massive 
halos with DM masses of 10 13-14 M , and are usually the cen - 
tral galaxies of rich groups and clusters (Zhe ng et al.||2 009). 
These high masse s have been confirmed thro ugh weak lensing 
measurements by Mandelba um et al.| ( [2006| ). |Ho et al.| |2009 ) 
have shown that X-ray luminous clusters with virial masses 
of 1 - 5 x 10 14 M at z ~ 0.2-0.6 contain LRGs, and that the 
brightest LRGs lie at the centers of the clusters. Other LRGs 
within a cluster, however, need not be found near a cluster's 
center. 

A simple explanation for the lack of Mg II from galaxies in 
our sample is that the gas in the halos of LRGs is too hot to 
sustain clouds cool enough to produce Mg II absorption. High 
mass halos are traditionally expected to gain gas through "hot 
mode" accretion, whereby gas from the intergalactic medium 
(IGM) is accreted into DM halos and is shock heated to its 
virial temperature, ~ 10 7-8 K for the most massive cluster ha- 



los (e.g.[Rees & Ostrikerj 19771 |White & Rees| 1978 [| White & 



r|1977[ 

Frenk 1991). Additional feedback mechanisms, such as su 
pernovae driven winds, or AGN accretion, can help maintain 
high gas temperatu res around ga laxies in massive halos (see, 
e.g. the review by Baugh 2006). For gas that is collision- 
ally ionized and in thermal equ ilibrium, most magnesium i s 
in the form of Mg III at 10 5 K ( [Sutherland & Dopita|[T993] ). 
Lower mass halos hosting ~ L* galaxies also have virial tem- 
perat ures that are too high for M g II to dominate the ion frac- 
tion (Fukugita & Peebles 2006), but models in which cool 
photoionized gas can survive in the hot halo are able to re- 
produce many of the characteristics of M g II (and other ion) 
systems (e.g. Mo & Miralda-Escude 1996). Semi-analytical 
models and numerical simulations of gas hydrodynamics in 
CDM constructs also find that cool gas can exist in galaxy 
halos through "cold mode" ac cretion, which occu rs when gas 
T < 2.5 x 10 5 K in|Keres et al.|d2009b] is funneled 



[with, e.g.. 

from the IGM down filaments into the halos. This pro cess is 
expected to dominate in low mass (< lO n ~ 12 M ) halos (|Birn 



9}]i 
rhis 



boim & Dekel 200 3"l|Katz et al.|2003l|Keres et al.|2005[|Eirn- 
boim et al. 12007 [IKeres et al.|2009b. 



Several investigators have considered whether the masses 
of Mg II absorbing halos are consistent with t hose expected 
for the existence of Mg II absorbing clouds. |Bouche et al. 
(2006j) cross-correlated Mg II absorption systems detected in 
SDSS spectra with a catalog of LRGs (with photo-zs), from 
which they were able to measure a mean halo mass 20 - 40 
times smaller than LRG halos. They found that the halo mass 
was anti- correlated with Mg II EW, and hence that the ab- 
sorbing clouds were not virialized in the halo. This result was 
unexpected, since more massive virialized halos, with brighter 
galaxies, should have absorbing gas clouds with larger veloc- 
ity dispersions. Since the EW of a Mg II line is, to first order, 
proportional to the number of narrow components comprising 
the line (at least for strong lines), then (assuming that Mg II 
absorption is directly related to clouds within a galactic halo), 



W r should be proportional to the virial velocity of gas in the 
halo, and hence to halo mass. Bouche et al. estimated the 
mass of halos in which strong (W r > 1 A) Mg II was found 
to be « 10 12 M , and determined that an anti-correlat i on ex- 
isted between halo mass Mh and W r . Lundgr en et al.| (f2009 ) 
refined these values with additional SDSS data releases, and 
determined that Mg II absorbers with W r > 1 .4 A arose in halo 
masses of 1O 1L3 M , while absorbers with 0.8 < W r < 1.4 A 
arose in halos with M h ~ 1O 127 M . 

Our results are in broad agreement with this anti-correlation 
of W r and Mh, in the sense that we find few Mg II systems in 
the high mass halos of LRGs. Bouche et al. did not predict 
a covering factor for a given halo mass, making it difficult 
for us to compare in detail our results with their models. In 
interpreting their results, Bouche et al. favored the explanation 
that galactic outflows within the halos caused the majority of 
Mg II systems, and not inf ailing IGM gas. In keeping with this 
hypothesis, we would say that we see little Mg II absorption 
because few star-forming galaxies exist in high mass groups 
and clusters. The LRGs in our sample are likely to be cluster 
ellipticals, with little significant star formation or outflowing 

c ool gas. We return to th is p oint below. 

Tinker & Chen] d2008] ) and |Chen & Tinker| ( |2008] ) modeled 



Mg II absorbing halos in terms of the HOD formalism, and 
found that in order to reproduce the anti-correlation of W r and 
Mh seen by Bouche et al., their models needed to incorpo- 
rate a transitional halo mass at which shock heated gas be- 
comes more dominant in higher mass halos, leading to less 
Mg II absorption. In these models, the halos of Mg II absorb- 
ing galaxies have masses of ~ 1O 11_13 M , with a transitional 
mass of 1O 1L6 M , a result consistent with the simulations of 
cold-mode accretion of gas into galaxy centers. Again, in the 
context of our survey, Mg II is not detected near LRGs be- 
cause the galaxies reside in halos that are too massive to sup- 
port cool Mg II clouds. 

In detail, it is less clear whether our results are compatible 
with Tinker & Chen's models, which were calculated for halo 
masses as high as 1O 14 M . They found that "a significant 
fraction of cluster- sized halos contain some cold gas with a 
high covering factor", although they noted that the constraints 
on their models were poor for halo masses > 1O 13 M . If we 
take all the LRGs in our sample to represent cluster-sized ha- 
los, our results would tend to argue against these findings. We 
would say that our sample shows that a significant fraction 
of cluster-sized halos contain cold gas with a low covering 
factor, or, alternatively, only a few cluster-sized halos contain 
cold gas with a high covering factor. 

4.3. Possible Origins for Observed Mg II Absorption 

Despite the fact that most LRGs do not cause Mg II absorp- 
tion, we have shown that even after removing a contribution 
to the number of detections that arise from chance along any 
sightline, there are still a small number of Mg II systems asso- 
ciated with the identified galaxies. In this section we discuss 
possible origins for these systems. 

Could the few Mg II systems we see around LRGs arise 
from the small amounts of cool gas that remain in galactic 
halos? Cold-mode accretion is expected to be a dominate 
channel for galaxies at high redshifts, or low-mass galaxies at 
lower redshifts. As indicated above, the halo masses of LRGs 
are likely to be several dex above the transition mass between 
the hot- and cold-mode regimes (1O 11_12 M ), so cold accre- 
tion is unlikely to be significant. At z = 0, numerical Simula- 
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tions suggest that the fraction of gas in halos with masses of 
1O 13 M and temperatures < 2.5 x 10 5 K is only a few per- 
cent, although whether cold filamentary streams reaching into 
the halo cor es at z ~ 0.5 are able to survive is unclear ( Keres 
|et al.|[2009] ). The problem with this speculation is that we 
do not actually know the mass distribution of the LRGs we 
have sampled. Without a better measurement of the LRGs' 
masses and environments, it is difficult to compare our results 
with simulations and determine whether even small amounts 
of cold gas might be associated with the LRGs we have stud- 
ied. 

A different explanation is that Mg II absorption might arise 
in the disks of interloping late-type galaxies which in habit the 
outer region s of the LRG groups and clusters (e.g. Hansen 
|et al.||2009] and refs. therein). Any difference in the red- 
shifts between an LRG and another cluster or group member 
will be much smaller than the error in the photometric red- 
shift of the LRG, and will not be detectable in Figure [9] The 
morphologies of galaxies are known to change from fate- to 
early-type as the distance to the center of a cluster decreases, 
and late-types may only need to present a small cross-section 
to account for the low covering factor seen in Figure [6] More- 
over, the transformations that take place in late-type galaxies 
in the environment of rich groups and clusters offer plausible 
mechanisms for increasing their Mg II cross-section; tidal in- 
teractions with other cluster galaxies, or stripping of galactic 
ISMs by the intracluster medium, can re-distribute cool gas 
over larger areas than those covered by single galaxy disks. 
The problem with this scenario, however, is that we would 
expect the surface density of late-type galaxies (as a function 
of the projected distance from the center of a cluster) to be rel- 
atively constant over Mpc scale^J there is no obvious reason 
why the number of spiral galaxies would decline to zero at im- 
pact parameters of 200 kpc from an LRG. Put another way, if 
late-type galaxies at the outer regions of galaxy clusters were 
responsible for detected Mg II systems, we would expect the 
number of systems "accidentally" intercepted by QSO sight- 
lines through a cluster to be largely constant between 20 and 
200 kpc of an LRG, possibly at a level slightly higher than the 
background Mg II systems discussed in ^3] Since this is not 
what we observe, we consider this scenario to be an unlikely 
explanation for the few Mg II systems we do detect. 

A more likely explanation is that the Mg II systems arise 
as a consequence of wet mergers between LRGs and gas-rich 
satellites, with the latter providing fuel for star formation in 
the former. If hierarchical merger models are correct, LRGs 
should certainly have the most vigorous merger histories of 
galaxies in the universe. Again, galaxy-galaxy interactions, 
or galaxy cannibalism, is a good way to distribute tidal debris 
over areas of the sky larger than those of individual galac- 
tic disks. Absorption from tidally disrupted gas would also 
explain the large Mg II EWs, since the kinematics of the in- 
teracting galaxies coul d spread Mg II com ponents over large 
velocity intervals (e.g. Bowe n et al.|1 994). More importantly, 
however, if Mg II systems are more closely associated with 
star-forming galaxies, as discussed in §1, where galactic out- 
flows push cool gas out of galactic disks, then perhaps LRGs 
that have undergone recent bursts of star formation may be 
responsible for strong Mg II absorption lines. Again, out- 
flows from star- forming galaxy disks w ould explain the strong 
Mg II EWs (e.g. Chelouche & B owen|20T0| and refs. therein). 



The star forma tion history of LRGs h as been studied by 
several authors. Eisen stein et al.| ( |2003| ) found that 10% of 
the S PSS LRG popula t ion ha ve significant [O II] emission, 
while IRoseboom et al.| 42006] ) found that 12.5% of L > 3L* 
LRGs have spectra with some combination of [O II] emission 
lines and strong US absorption lines, indicative of ongoing, 
or recent, star formation. The similarity between this frac- 
tion of LRGs, and the fraction we find associated with Mg II 
absorption, is striking. Excluding their [O II] emission-only 
subset of galaxies, ~ 4% of Roseboom et al.'s LRGs show 
good evidence for star formation; as noted by both Eisen- 
stein et al. and Roseboom et al., the detection of emission 
lines alone without corroborating absorption features is am- 
biguous, since [O II] can arise from galaxies with AGN activ- 
ity (such as LINERS), rather than from star formation. The 
Mg II absorption cross-section of AGN has not been stud- 
ied in detail, although QSOs at similar redshifts to the LRGs 
in our sample are known to show high Mg II covering fac- 
tors (~ 90%) out to ~ 1 00 kpc for lines with W r > 0.5 A 
( |Bowen et al.]2005||2mi) . Of course, QSO and AG N activity 
may itself b e triggered by galaxy mergers (see, e.g. Hopkins 
et al. 12008] and refs therein), a process which distributes gas 



over a wide cross-secti on, triggers high star formation (Mi 
hos & He rnquist 1996), and causes gas to be expelled from 



the ga laxy through supernovae driven winds ( Springel et al. 



2005). Hence, we may be seeing Mg II from the same com 
bination of star-forming and (perhaps) AGN LRGs as those 
selected by Roseboom et al. and Eisenstein et al., leading to 
the same fraction of Mg II absorbers. Roseboom et al. were 
unable to settle on the physical mechanisms which might pro- 
duce star-forming/AGN LRGs, but noted that their existence 
was consistent with ACDM cosmologies that included merg- 
ers of gas-rich galaxies, and cooling at the centers of clusters, 
much as already discussed above. Obviously, whether the 
LRGs in our sample show any signs of the same star formation 
or AGN activity seen by Eisenstein et al. and Roseboom et al. 
can be tested by sim ply recording the spectra of the LRGs. 
Gauthier et al. ( 2010) have begun a spectroscopic follow-up of 
their LRGs, but so far, only one galaxy associated with Mg II 
absorption and at an impact parameter < 100 kpc has been 
observed. Its SED shows no obvious signs of star-formation 
or AGN activity. 

Whether Mg II absorption originates in tidal debris from 
merging galaxies, or whether it arises in outflows from star 
formation instigated by merging galaxies, in both cases the 
processes involve the interaction of galaxies. Unfortunately, 
at redshifts of ~ 0.5, measuring the rate of even major LRG 
mergers, where both interactors might be bright enough to be 
visible in large-area sky surveys, has proved difficult; the rate 
though, does appear to be low ([ Masjedi et al. 2008; Bundy 
et al.|2009[|De Propris et al.|2010[ and rets therein). The co- 
moving volume density of the most massive (M >~ 1O U M0) 
galaxies changes little at z ~ 0.2- 1. The numbers a re not 
increasing [at least to within factors of 2-3 — see, e.g. Con- 
|selice et aL] ( |2007| )], even though more would be expected to 
arise at lower-z if they continued to accrete galaxies through 
mergers. 

At much lower redshifts, z — 0.02-0.14, where gal axy en- 



vironment is easier to map in shallow galaxy surveys, Huang 
&~Gu] ( [2009] ) found that 5% of bright (between M r of -20 



8 For example, the mean radius at which the galaxy density is defined to 
be 200 times the field density CR200) is usually given to be ~ 1 Mpc. 



and -24, so similar to the LRGs) E and SO galaxies showed 
signs of ongoing or recent star formation, and an additional 
10-13% showed LINER or Seyfert activity. These numbers 
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are similar to those measured specifically for the LRGs by 
Roseboom et al. Unfortunately, Huang & Gu excluded galax- 
ies with disturbed morphologies, making it hard to know how 
much star formation in early-type galaxies might be associ- 
ated with merging systems. Interestingly, the environments 
of their star-forming ellipticals appear to be quite mixed. Al- 
though 5 galaxies seem to be relatively isolated, 8 of the 13 
listed by Huang & Gu have one or more neighbors with simi- 
lar redshifts within 300 kpc and/or have been associated with 
galaxy groups or clusters by other authors. Hence, at low-z 
it seems plausible that a small number of bright, star-forming 
ellipticals can be found in similar dense environments to those 
of LRGs at higher-z. 

Can we measure any difference in the environment of those 
LRGs that we associate absorption with and those that show 
no absorption? To investigate this question, we collated all 
galaxies within 5 arcmin of every QSO sightline in our survey 
from the SDSS DR7 GALAXY database. At z = 0.6, 5 arcmins 
corresponds to a proper length of 2 Mpc, which approximates 
cluster-like scales. We calculated the number of galaxies per 
shell along each sightline, n(r), where we counted galaxies 
in shells of width 50 arcsec, and then compared both the av- 
erage and median of n(r) for all sightlines with and without 
Mg II detections. We took the error in these values to be the 
measured standard deviation in n(r). In one case we counted 
all galaxies, regardless of their redshifts; in a second cas e, we 
used only red g alaxies that satis fied the color selection of Can- 
non et al.| ( |2006| ) discussed in j j2.3| and shown in Figure [3] In 
the latter case, these red galaxies are more likely to be at the 
redshifts of the previously selected LRGs. 

For both galaxy samples, unfortunately, we could not differ- 
entiate between the average or median of n(r) for sightlines 
with and without Mg II absorption. We found that the vari- 
ance in the number of galaxies per field dominates the distri- 
bution of n(r) for both absorbers and non- absorbers. In par- 
ticular, having only 85 unique QSO sightlines showing Mg II 
detection provides too small a sample to show any difference 
with the distribution of galaxies in the 482 QSO fields that 
show no Mg II. A larger study which matches all LRG envi- 
ronments around all QSO sightlines, and not just QSOs with 
LRGs within 30 arcsec of a sightline (as used in this study), 
may be more fruitful, but is beyond the scope of this paper. 

Our search for Mg II absorption lines from LRGs can only 
be used so far in interpreting the origin of QSO absorption 
line systems. The lack of Mg II from this population of galax- 
ies suggests that the majority of Mg II systems have little to 
do with massive red galaxies, nor probably, massive hot halos. 
The implication is that if Mg II systems are indeed associated 
with galaxies (and are not simply metal-enriched intergalactic 
clouds) then other types of galaxies and environments are pri- 
marily responsible for Mg II absorption. This is hardly a sur- 
prising conclusion, given the identification of the strong links 
between Mg II systems and late-type, star-forming galaxies 
discussed in the Introduction of this paper. Indeed, as our un- 
derstanding of galaxy properties and galaxy number counts 
has become more refined, and with a more precise counting 
of Mg II systems at the lowest redshifts, it has become easier 
to account for the origi n of all Mg II systems by s tar form- 
ing disk galaxies (e.g. Chelouch e & Bow en 2010, and refs 
therein). Whether the small percentage of Mg II absorbers 
found herein to be associated with LRGs provides any deeper 
insights into the origin of Mg II systems, is not clear. 
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